Joule-Thomson coefficient of ideal anyons within fractional exclusion statistics 
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The analytical expressions of the Joule-Thomson coefficient for homogeneous and harmonically 
trapped three-dimensional ideal anyons which obey Haldane fractional exclusion statistics are de- 
rived. For an ideal Fermi gas, the Joule-Thomson coefficient is negative, which means that there is 
no maximum Joule- Thomson inversion temperature. "With careful study, it is found that there exists 
a Joule- Thomson inversion temperature in the fractional exclusion statistics model. Furthermore, 
the relations between the Joule-Thomson inversion temperature and the statistical parameter g are 
investigated. 
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I. INTRODUCTION 

The fractional exclusion statistics of Haldane is an in- 
termediate statistics between Bose-Einstein and Fermi- 
Dirac statistics. It is also a generalized dimensional- 
independent statistics based on state counting meth- 
ods and is suitable to describe interacting many-particle 
systems in condensed matter [H-Q. By considering TV 
particles in a c?-dimensional Hilbert space at fixed size 
and boundary conditions, the linear relation between the 
changes of the single-particle space Ad and the changes 
of the particle number AA^ is defined as Ad = —gAN 
with a parameter g given by Haldane [ij. 

In Haldane fractional exclusion statistics, the number 
of microscopic quantum states W of N identical particles 
occupying a group of G states is 



w = Y[ 



[Q + (jV.-l)(l-g)]! 
m[G, - gN, - (1 - gW 



(1) 



which is interpolated by Johnson, Canright and Wu 
It corresponds to the result of Bose-Einstein statistics 
when the weight factor g — and to the one of Fermi- 
Dirac statistics when <? = 1. Here the statistical parame- 
ter g is described as the change in the number of available 
states when one particle is added to the system. 

The anyon statistical model has been extensively stud- 
ied in the literature [3-[TT| . The thermodynamic solution 
of the one-dimensional ideal anyon gas which obeys the 
Haldane statistics is equivalent to the Bethe ansatz so- 
lution of the Calogero-Sutherland model [T2| - [l5j . Fur- 
thermore, the thermodynamic extension to d ^ 1 di- 
mensions of the Calogero-Sutherland model is explored 
by Potter et al. (3). The Haldane's fractional exclusion 
statistics is widely used to describe the low-dimensional 



condensed matter physics. In the literature, the prop- 
erties of spinous are characterized in terms of the Hal- 
dane's fractional exclusion statistics with corresponding 
models ( e.g., the Haldane- Shastry model in which free 
spinous exist [lB-[l3 , the one-dimensional supersymmet- 
ric t- J model in which the thermodynamics of spinous is 
investigated and the Wess-Zumino-Witten model 

Besides, the Hubbard model with an infinite- 
range interaction is also used to study particles obeying 
Haldane fractional statistics [1^ [2^ . 

In thermodynamics, the fact that the temperature 
changes with the decrease of pressure during an adia- 
batic or throttling expansion is called the Joule-Thomson 
effect. The adiabatic expansion or Joule- Thomson pro- 
cess describes the procedure that a gas is forced through 
a porous plug without heat exchange with the environ- 
ment. This effect can be described by the Joule-Thomson 
coefficient ujt, which is the partial derivative of temper- 
ature with respect to pressure at constant enthalpy 



UjT 



dT 
dP 



(2) 
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where P is pressure, T is system temperature and H is 
enthalpy. If ujt > 0, it shows that the temperature 
of the system decreases during the adiabatic expansion. 
On the other hand, if ujt < 0, the temperature increases 
in the throttling process. By setting ujt = 0, one can 
calculate the maximum Joule-Thomson inversion tem- 
perature. Below this inversion temperature, a free adia- 
batic expansion causes a decrease in temperature, while 
it causes a temperature increase above this inversion tem- 
perature [27| . 

There have been many works on the Joule-Thomson 
coefficient. The Joule- Thomson coefficients of ideal quan- 
tum systems are obtained by Ref. [28]. The corre- 
sponding results for weakly interacting Fermi and Bose 
gas es are given by means of the pseudopotential method 
[291 . Recently, the Joule-Thomson coefficient for a d- 
dimensional ideal Bose gas was derived in a power-law po- 
tential [30| . The Joule- Thomson coefficient for a strongly 



2 



interacting Fermi gas was also discussed within the quasi- 
hnear approximation framework j3T| . 

According to Landau's phenomenological theory, there 
is an important parameter called effective mass which is 
used to characterize the quasiparticle and is determined 
by the low-temperature isochore heat capacity. If ideal 
anyons are considered as quasiparticles to describe inter- 
actions j32h36i] , the effective mass is a significant physical 
quantity in the fractional exclusion statistics model. 

In this paper, the analytical expressions of Joule- 
Thomson coefficients for homogeneous and harmoni- 
cally trapped three-dimensional ideal anyons are de- 
rived within the Haldane fractional exclusion statistics. 
Furthermore, the Joule-Thomson inversion temperature 
varying with the statistical parameter g is plotted and 
the effective mass is calculated. The outline is as fol- 
lows. In Section ini the distribution function of Haldane 
fractional exclusion statistics is obtained by the method 
of Lagrange multiplier. The Joule-Thomson coefficients 
of homogeneous and harmonically trapped ideal anyons 
within fractional exclusion statistics are derived analyti- 
cally in Section Um Effective mass is calculated in Section 
IIVI The numerical calculations and the Joule-Thomson 
inversion temperature plotted as a function of parameter 
g are obtained in Section [V] A summary is given in the 
final section. 



II. DISTRIBUTION FUNCTION OF HALDANE 
FRACTIONAL EXCLUSION STATISTICS 

Let us set two Lagrange multipliers a = —fJ./T and (3 = 
l/T, where fi is the chemical potential and the natural 
units ks = h = 1 are used. The average occupation 
number is defined as Ni = Ni/d. According to the 
Lagrange multiplier method 6\nW — aSN — /3SE = 0, 
the most probable distribution function can be derived 
as (3 



N 



1 



(3) 



where the statistical parameter g and uj are related as 



e = M + T[(l-5)ln(l 



g In cj] 



(4) 



where e is the single-particle energy. When g = 1, the dis- 
tribution function N goes back to fermionic distribution 
function, and when g = it becomes a bosonic one. Ac- 
cording to equation ([3]), when the temperature is zero, 
one can get N = with e > fi, and N = 1/g with 
e < ^, which indicates that g characterizes the gener- 
alized Pauli principle since the maximum value of the 
occupation number for a single-particle state is l/g. 

We define the fugacity z in terms of the chemical po- 
tential fi and introduce a parameter ujq as follows: 



exp 



'-'ton' 



(5) 



Inserting equation ([5]) into equation ([4]), one obtains 

1 +CJ 



T 



(1-5) In 



gin 



(6) 



III. THE JOULE-THOMSON COEFFICIENT 
GIVEN BY THE FRACTIONAL EXCLUSION 
STATISTICS 

A. Homogeneous gas 

The density of states is D{e) = {2mf/'^Ve^^^/{2n'^) 
with two degrees of the spin degeneracy for a homoge- 
neous gas, where m is the particle mass and V is the 
system volume. 

The_ expression for the grand thermodynamic potential 
il is 



PV 



TV{2m) 
3^ 



3/2 



^3/2 In ( 1 + ± 

UJ 



(7) 



In the one-dimensional supersymmetric t-J model, three 
species of particles will contribute to the grand thermo- 
dynamic potential by applying an external magnetic field, 
respectively [20j. In the literature, it is pointed out that 
the one-dimensional supersymmetric t-J model is equiv- 
alent to Haldane fractional statistics with properly fixed 
statistical parameters [l^. Here, we limit to the three- 
dimensional free particles in terms of the Haldane's frac- 
tional exclusion statistics without magnetic field. The 
corresponding thermodynamic potential is given by equa- 
tion ©. 

With the help of equation ([6]), equation ([7]) can be 
reduced to 



2T 

P = -jjG5/2{z,g), 



(8) 



where the thermal de Broglie wavelength is defined as 
A = ^2TT/{mT) and 



G„{z,g) = 



r{n) 



9 



(9) 



is the Calogero-Sutherland integral function f3|. T{n) = 
{n — 1)! = Jp°° exp{—y)y^^^^dy is the gamma function and 
X = e/T. The Calogero-Sutherland integral function has 
been proved to satisfy the recurrence relation which is 
the same as the Bose-Einstein and Fermi-Dirac integral 
functions Ijl. The recurrence relation is 



d 

z — Gn{z,g) = G„_i(z, ( 



(10) 
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By turning the variable e into uj through equation ([6]) , 
equation ^ can be reduced to 



Gn{z,g) 



r(n 



(11) 



where 



hn{ujo,g) 



In 
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UJ{UJ + 1) 



(12) 



The finite-temperature particle number density and in- 
ternal energy density can be represented as 



E 
V 



l_ r°° D{e)de 
V Jo ^ + 9 
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(13) 



uj + g 



-G5/2{z,9) 



^ -P. 



(14) 



At zero temperature, the particle number is = 



il/g)J^^ D{e)de = i2mEF)^/^V/{3T:^), where Ep is 
the uniform ideal Fermi energy and Ep satisfies Ep — 
g^/^Ep. By replacing the particle number of the ground 
state into equation (fT3|) . one can get 



3^1/2 / T 



3/2 



03/2(^,5) = 



T 
2 VT> 
1 



3/2 



hi/2{^o,9) 
(15) 



with the Fermi characteristic temperature Tp for a uni- 
form ideal Fermi gas. 

From equation (jl3p . the partial derivative of particle 
number density n with respect to temperature T at con- 
stant n is written as 



dn 
df 



( 5G3/2(z,g) 



7^3/2(2:, 5) 



T 



T 



V dT 



(16) 



By combining equation ((T^ with {dn/dT)n ~ 0, the 
partial derivative of chemical potential to temperature 
takes the form 



dji\ _^ 3^3/2(2:,. 9) 
dTj^^''' 2Gi/2(z,.9)- 



(17) 



According to the constant total particle number N ^ 
one can have 



dP 
df 



dV 

df 



dP 
df 



1 

A3 



5G5/2(z,5) 



3G^/2(^>g) 
^1/2(^,5) 



(18) 



\ dT 



ATG3/2{z,g) 

^Gz/2{z,g)Gi/2{z,g) 
^3/2(^75) 



(19) 



To derive the expression of the Joule- Thomson coef- 
ficient, the isochore heat capacity Cy and isobar heat 
capacity Cp per particle are derived first as 



N 



1 /dE 

N \dfyj,y 

15^5/2(2, g) _ 9^3/2(2, ff) 
4G3/2(z,5) 4Gi/2(z,5) 
5/13/2(^0,5) 9/ii/2(a;o,5) 



Gp 
N 
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cv T /ap\ /ay' 
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5/i3/2(^0,g) 

6/11/2(^0,5) 
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(20) 
(21) 
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(22) 



Furthermore, from the fundamental thermodynamic 
relations and equations ©, and (f^ . the Joule- 

Thomson coefficient can be given by 



UJT 



T 



1 



dV_ 
'df 
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(23) 



B. Harmonically trapped gas 



We define the geometric mean of the trap frequencies 
as tn = (wxtOytOzY^'^ ■ The corresponding density of states 
is D{e) = /w^ . In the similar way as described in the 
homogeneous case, the Joule- Thomson coefficient of a 
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trapped gas is given as 
V /n7\3 



UJT 



6" \t) 



1 



Giiz,g) 
G2{z,g) 



4Gi{z,g)G2iz,g)-SGl{z,g) 



1 

r3 



2hi{u}o,g) 



8h3{uJo,g)hi{u;o,g) - 9/i|(wo,5) 



(24) 



with V = glS^. 

The corresponding analytical expressions for the par- 
ticle number and isobar heat capacity per particle are 



iV = 2(^) G3{z,g), 



Cp _ AhsjujQ^g) 
N ' 3/^2^,5) 



8h3{uJo,g)hi{uJo,g) 
3/il(wo,5) 



At zero temperature, the particle number is 



TV = 



1 fEf 



(25) 



(26) 



(27) 



3 \ zu 

Substituting equation (|27)) into equation ([25)) . one gets 

= 1. (28) 



IV. EFFECTIVE MASS 

At very low temperature, the G„ (z, g) can be expanded 
as I, [Mil 



Gn{z,g) 



(Inz)" 
5r(n+l) 



1 



TT^ (7n(n — 1) 
"6 (lnz)2 



(29) 



From equations (PO]) and (P^ . the expression of the 
isochore heat capacity at extremely low temperature can 
be given by 



N 2 ^ VT'f 



The zero-temperature effective mass is 

Cv 



m 
m 



{Cv)ide 



9 



1/3 



where {Cv)ideai ~ tt'^TN/ {2Tp) is the isochore heat ca- 
pacity of an ideal Fermi gas at the first-order approxima- 
tion and m is the mass of non-interacting ideal fermions. 



V. DISCUSSIONS 

The numerical results will be given in this section. 

The Joule- Thomson coefficient plotted as a function 
of reduced temperature can be obtained from equations 
(fist and (l23t for a homogeneous gas of anyons. It can 
also be given by equations (|24|) and ([28|) for a trapped 
anyon gas numerically. 
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FIG. 1: The Joule-Thomson coefficient ujt is plotted as a 
function of the reduced temperature for homogeneous gases. 
The solid curve denotes that for the ideal fermions, and the 
dashed curves denote the ones for ideal anyons with different 
statistical parameter g = 0.1, ^ , 0.5, 0.7, 0.9. Inset shows that 
the Ujt with fixed g = -^f changes its sign at T/Tp ~ 0.4695. 
m — 1 and Ep = 1 are chosen for convenience. 




FIG. 2: The line styles are similar to the ones in figure [T] for 



harmonically trapped gases. The ujp with g 



changes 



(30) its sign at T/Tp x 0.3369 in the inset. Here, Ep = 1 

Figures [T] and [H show that the Joule- Thomson coeffi- 
cient is larger for smaller values of g at the same tem- 
perature for both homogeneous and trapped anyons. Be- 
sides, the curves all tend to zero in the high temper- 

(31) ature limit as a Boltzmann gas and approach the mi- 
nus infinity in the strong degenerate limit T — ^ 0. The 
Joule- Thomson coefficients for both homogeneous and 
trapped ideal Fermi gases are always negative, which 
means that there is no Joule- Thomson inversion tempera- 
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ture for ideal Fermi gas. During the adiabatic expansion, 
the temperature always increases for an ideal Fermi gas. 
For ideal anyons, however, there exists a Joule-Thomson 
inversion temperature, which depends on the parameter 
g in the fractional exclusion statistics model. 

Recently, the strong interaction in two-component ul- 
tracold fermions is a hot topic [38l - l43| . The strongly in- 
teracting Fermi gas is called the unitary Fermi gas [iol - 
|43| . As a hypothesis, the three-dimensional ideal anyons 
with fractional exclusion statistics can be used to model 
the statistical behavior of a real unitary Fermi system 
[s^ - lscj . This fractional exclusion statistics hypothesis 
is found to be in good agreement with the experimental 
results in a harmonic trap [sBl- 

The parameter q in the statistical model can be fixed 
as g = ^^/^ = ^ |35l . [36| , where the universal constant 
^ = I is given by the developed quasi-linear approxi- 
mation j43 - l49| . By taking g = ^ as an example, the 
inversion temperature of a homogeneous unitary Fermi 
gas is T/Tp w 0.4695, and it is T/Tp « 0.3369 for a har- 
monically trapped unitary system as indicated in figures 
□ andm 
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FIG. 3: The reduced Joule-Thomson inversion temperature 
T/Tf versus g. The dashed curve denotes the result of ho- 
mogeneous anyon gas, and the solid curve represents that of 
trapped gas. 

It is evident from figure [3] that the reduced Joule- 
Thomson inversion temperature increases with the in- 
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